The ortho-to-para ratio of ammonia in the LI 157 outflow 
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ABSTRACT 

We have measured the ortho-to-para ratio of ammonia in the blueshifted gas 
of the L1157 outflow by observing the six metastable inversion lines from (J, K) 
— (1, 1) to (6, 6). The highly excited (5, 5) and (6, 6) lines were first detected 
in the low- mass star forming regions. The rotational temperature derived from 
the ratio of four transition lines from (3, 3) to (6, 6) is 130-140 K, suggesting 
that the blueshifted gas is heated by a factor of ~10 as compared to the 
quiescent gas. The ortho-to-para ratio of the NH3 molecules in the blueshifted 
gas is estimated to be 1.3-1.7, which is higher than the statistical equilibrium 
value. This ratio provides us with evidence that the NH3 molecules have been 
evaporated from dust grains with the formation temperature between 18 and 25 
K. It is most likely that the NH3 molecules on dust grains have been released 
into the gas phase through the passage of strong shock waves produced by the 
outflow. Such a scenario is supported by the fact that the ammonia abundance 
in the blueshifted gas is enhanced by a factor of ~5 with respect to the dense 
quiescent gas. 

Subject headings: ISM: individual (LI 157) — ISM: jets and outflows — ISM: 
molecules — stars: formation — shock waves 
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1. INTRODUCTION 

The inversion lines of metastable ammonia have been extensively used for observations 
of dense (~10^ cm~^) molecular cloud cores (e.g. Myers & Benson 1983; Benson & Myers 
1989), while the fractional abundance of NH3 varies among the clouds (Benson & Myers 
1983). On the basis of a systematic survey, Suzuki et al. (1992) pointed out that the NH3 
tends to be more abundant in the star-forming cores than the starless cores. Such a trend 
is interpreted in terms of the gas-phase chemical evolution which predicts that the NH3 is 
deficient in the early stages of chemical evolution and becomes abundant in the later stages 
(e.g. Herbst & Leung 1989; Suzuki et al. 1992). In addition to the gas-phase chemical 
evolution, the gas-grain interaction is also considered to be important for the abundance 
variation of NH3 (e.g. d'Hendecourt, AUamandola, & Greenberg 1985; Brown & Charnley 
1990; Nejad, WiUiams, & Charnley 1990; Flower, Pineau des Forets & Walmsley 1995). In 
the star-forming dense cores, the NH3 molecules retained on dust grains can be released into 
gas phase by means of shock waves caused by the outflow and the radiation from the newly 
formed star (e.g. Nejad et al. 1990; Flower et al. 1995), and are expected to contribute to 
increase the NH3 abundance. 

Usually, it is not easy to estimate contribution of molecules desorbed from dust grains 
from those formed in the gas phase. However in the case of NH3, an ortho-to-para ratio of 
NH3 provides us with such information. Ammonia has two distinct species, the ortho-NH3 
{K = 3n) and the para-NH3 {K ^ 3n), which arise from different relative orientations of the 
three hydrogen spins. The ortho-to-para ratio is expected to be the statistical equilibrium 
value of 1.0 when the NH3 molecules are formed in the processes of gas- phase reactions 
or grain surface reactions. This is because the reactions to form the NH3 molecule release 
large excess energies compared to the energy difference between the lowest ortho and para 
states, which make the ortho-to-para ratio close to the statistical value. However, the 
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ortho-to-para ratio is expected to be larger than unity when the NH3 molecules adsorbed 
on grain surface released into the gas phase with the excess energy for desorption which 
is comparable to the energy difference between the ortho and para states. This is because 
the lowest energy level of the para species is 23 K higher than that of the ortho species, 
para species require more energy for desorption than ortho species. The time scale of the 
interconversion between ortho-NHa and para- NH3 is considered to be the order of 10® yr 
in the gas phase (Cheung et al. 1969). Therefore, the ortho-to-para ratio provides us with 
valuable information on physical conditions and chemical processes when the NH3 molecules 
are released into the gas phase. 

Since the lowest energy level of ortho-NHa, (J, K) — (0, 0), has no inversion doubhng, 
it is necessary to observe the transitions higher than {J, K) = (3, 3) to measure the 
ortho-to-para ratio. Such high transitions are hardly excited in the dark clouds of low 
temperature, however, recent observations have revealed that shocked gas associated with 
molecular outflows in the dark clouds is heated enough to excite the transitions higher 
than {J,K) = (3, 3). One of the prototypical objects with shock-heated gas is the bipolar 
outflow in the dark cloud LI 157 (e.g. Umemoto et al. 1992) which is located at 440 pc 
from the sun (Viotti 1969). Previous NH3 observations have revealed that the gas in the 
outflow is heated to more than 50-100 K (Bachiller, Martin-Pintado, & Fuente 1993; Tafalla 
& Bachiller 1995). A good morphological coincidence of the NH3 distribution (Tafalla & 
Bachiller 1995) with that of the SiO (J=2-l) emission (Zang et al. 1995; Gueth et al. 
1998), which is considered to be a good tracer of the shocked molecular gas, reveals that 
the hot ammonia arises from the shocked gas. The NH3 abundance enhancement observed 
in the shocked gas (Tafalla & Bachiller 1995) suggests the possibility that the NH3 retained 
on grain mantles is released into the gas-phase by means of shocks. 

In this Letter, we report the observations of six metastable inversion lines of NH3 from 
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{J,K) — (1, 1) to (6, 6) toward the blueshifted lobe of the LI 157 outflow. We have detected 
the high excitation NH3 (5, 5) and (6, 6) Unes for the first time in the low-mass star forming 
regions. Of six observed transitions, the (3, 3) and (6, 6) states belong to ortho- NH3 and 
the other four states belong to para-NHs. The detection of both (6, 6) and (3, 3) emission 
enables us to measure the ortho-to-para ratio which provides us with information on the 
contribution of ammonia molecules desorbed from grains. 

2. OBSERVATIONS 

The observations were carried out in 1991 June and 1992 May with the 45m telescope of 
the Nobeyama Radio Observatory. We observed six inversion transitions of the metastabel 
ammonia, {J,K) — (1,1), (2, 2), (3, 3), (4, 4), (5, 5), and (6, 6), toward three positions 
along the blueshifted lobe of the CO outflow. The position offsets from IRAS 20386+6751 
(hereafter referred to as L1157 IRS), a(1950) = 20^38'"39.'^6, 5(1950) = 67°51'33", are 
(0", 0"), (20", -60"), and (45", -105"), which were referred as positions A, B, and C, 
respectively by Mikami et al. (1992). Positions B and C correspond to the tips of two CO 
cavities (Gueth, Guilloteau, & Bachiller 1996), toward which strong NH3 (3, 3) and SiO 
( J=2-l) emission lines are observed (Tafalla & Bachiller 1995; Zhang et al. 1995; Gueth 
et al. 1998). We also observed the (1, 1), (2, 2), and (3, 3) transitions toward additional 
two positions outside of the lobe and one position toward the peak of the redshifted CO 
emission. The observed positions are shown in Figure 1 superposed on the CO (J=l-0) 
outflow map observed by Umemoto et al. (1992). The lowest three transitions were 
observed simultaneously. The higher three transitions were also observed simultaneously. 
At the observed frequencies, the telescope had a beam size of 72" and the beam efficiency 
of 0.8. The frontend was a cooled HEMT receiver whose typical system noise temperature 
was ~200 K. The backend was a bank of eight acousto-optical spectrometers with 37 kHz 
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resolution, which corresponds to a velocity resolution of 0.50 km s~^. All the observations 
were made in position switching mode. To obtain the spectra, we have integrated 45 
minutes per position, the resulting rms noise level per channel was ~ 0.02 K. 

3. RESULTS 

In Figure 2, we show the spectra of six transitions observed at positions A, B, and C. 
At position A, a strong (1, 1) line and weak (2, 2) and (3, 3) lines are observed. No emission 
from the transitions higher than (4, 4) has been detected. The (1, 1) line observed at 
position A shows the quiescent component peaked at Vlsr = 2.85 km s~^, which is almost 
the same as the cloud systemic velocity of ~2.7 km s~^ derived from the ^^CO observations 
(Umemoto et al. 1992), with a weak blueshifted wing. The (2, 2) and (3, 3) hues appear 
in the velocity range of the (1, 1) wing emission, suggesting that these line profiles are 
contaminated by the blueshifted component in the southern lobe because of the large beam 
size of 72". Since the (2, 2) and (3, 3) lines show no significant emission at the velocity of 
the quiescent component, the rotational temperature. Trot, of the quiescent component is 
estimated to be ^10 K. 

At positions B and C, all six transitions up to (6, 6) were detected. The (5, 5) and (6, 
6) lines are first detected in the low-mass star forming regions. All transition lines except 
the (1, 1) one have peak velocities blueshifted by 1-2 km s~^ from the cloud systemic 
velocity and show broad widths of AV ~4-10 km s"^ (measured at la level), indicating 
that the emission of the higher energy level arises from the high velocity gas. The (3, 3) 
line profiles at positions B and C resemble the profiles of the SiO (J=2-l) and CS (J=2-l) 
hues observed by Mikami et al. (1992). Among the six transition hues, the (3, 3) line is 
the strongest as previously pointed out by Bachiller et al. (1993). The peak brightness 
temperatures of the (2, 2), (3, 3), and (4, 4) lines observed at position B are lower by a 
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factor of 1.5-2 than those observed by Bachiller et al. (1993) at the same position. This 
may be due to the beam dilution effect, because the spatial extent of the ammonia emitting 
area is smaller than our beam ( Bachiller et al. 1993; Tafalla & Bachiller 1995). 

Toward three other positions outside the blueshifted lobe, we have detected only the 
(1, 1) emission line at the systemic velocity. 

4. DISCUSSION 

4.1. Temperature and Ortho-to-peira Ratio of the High- Velocity Gas 

To obtain the rotational temperature T^^t of the blueshifted gas, we assumed optically 
thin emission, and constructed rotation diagrams, i.e., the logarithms of NH3 column 
densities divided by statistical weight of the transition, plotted against the energy above the 
ground state. The contribution of the quiescent component in the (1,1) line was eliminated 
by performing the multi-components gaussian fitting. The rotational temperature estimated 
from the (1, 1) and (2, 2) data, Ti.ot(l,l; 2,2) is 38 K at positions B and is 48 K at position 
C. These temperatures are a factor of 4 to 5 higher than that of the quiescent component 
at position A. Figure 3 shows that the Trot obtained from the data of higher transitions 
are significantly higher than rrot(l,l; 2,2). The rotational temperature obtained from the 
para-NHg (4, 4) and (5, 5) data, Tj.ot(4,4; 5,5), are ~130 K for both positions B and C. The 
slopes between the ortho-NHs (3, 3) and (6, 6) data are almost parallel to those between 
the para-NHs (4, 4) and (5, 5) ones, indicating that Trot(3,3; 6,6) is comparable to Ti.ot(4,4; 
5,5). 

It should be noted that the column densities of the ortho species are higher than 
those of para species; this suggests that the ortho species are more abundant than the 
para species. If we assume an ortho-to-para ratio to be ~1.5, the data of four transitions 
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involving the ortho- and para-states align on the straight hnes; the best fits provide an 
ortho-to-para ratio at position B to be 1-7^° 3 and that at position C to be 1.3±0.2. Then 
we obtain the T^t at position B is 140]|]3K (the upper and lower limit of the T^ot correspond 
to those of the ortho-to-para ratio, respectively) and that at position C is 125 ['^^K. 

The discrepancy between the Trot(l,l; 2,2) and the Trot derived from the higher 
transition data can be explained by two components of gas with different kinetic 
temperatures as suggested by Avery & Chiao (1996) from their SiO observations. However, 
we consider that the most of the gas traced by the NH3 emission is heated to 130-140 K 
in the following reason. It is known that the metastable populations may deviate from a 
true Boltzmann distribution due to coUisional depopulation of the higher nonmetastable 
levels. As argued by Danby et al. (1988), the Trot tends to underestimate the gas kinetic 
temperatures T^, at the T^ range higher than ~ 30 K. Since the difference is remarkable for 
the Trot derived from the lower J transitions, the Trot estimated from the higher transitions 
is considered to be the better indicator of the kinetic temperature than Trot (1,1; 2,2). Such 
two rotational temperatures were also observed in the M17SW molecular cloud (Giisten & 
Fiebig 1988). 



4.2. Ammonia Abundance Enhancement 

The beam-averaged NH3 column density in the quiescent gas at position A is estimated 
to be A'"(NH3) = 2 x 10^^ cm~^. When we employed the H2 column density of the ambient 
gas derived from the ^'^CO and C^^O observations (A^(H2) = 3x10^^ cm^^), the NH3 
abundance for the quiescent component is estimated to be X(NH3) — 7 xlO~^, which is 
comparable to those measured in nearby dense cores (X(NH3) = (3-10) xlO"^; Benson & 
Myers 1983). The A^(NH3) in the blueshifted gas averaged over the beam at positions B 
and C are calculated to be 1x10^"^ cm~^ and 7x10^^ cm~^, respectively. The H2 column 
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densities of the blueshifted gas were estimated from the CO {J— 1-0) data obtained at the 
NRO 45m telescope (beam size was 16") by convolving the data with the 72" gaussian 
beam and assuming the optically thin CO emission, Tex = 130 K, and H2/CO abundance 
ratio of 10^. They are estimated to be 7x10^^ cm~^ toward position B and 2x10^^ cm~^ 
toward position C. By using these H2 column densities, we obtained ^(NHs) = lxlO~^ 
at position B and ^(NHa) = 3x10^^ at position C, which arc a factor of 2-5 higher than 
that of the quiescent component. The NH3 column densities derived by Tafalla & Bachiller 
(1995) from their high-resolution VLA data are ~5xl0^^ cm~^ toward both positions B 
and C. If we compare these values with the A'^(H2) in the 16" beam derived from the CO 
(J=l-0) data by assuming an excitation temperature of 130 K, which are 2x10^^ cm~^ for 
position B and 8x10^0 cm-^ for position C, we obtained ^(NHa) = 3x10"^ and ^(NHa) 
= 7x10"^ for positions B and C, respectively. The H2 column densities used here may be 
somewhat underestimated because the sizes of the NH3 enhanced regions (< 10") in the 
map of Tafalla & Bachiller (1995) are smaller than the beam size of the CO data. When 
we take this into account, the ammonia abundances obtained from the higher resolution 
data are consistent with those from the 72" resolution data. Therefore, we conclude that 
the NH3 abundance in the shocked regions is enhanced by a factor of ~5. 

4.3. Ortho-to-pcira Ratio and Its Implication for the Contribution of 

Desorbed Ammonia 

The derived ortho-to-para ratio in the blueshifted gas which is larger than the statistical 
value of 1.0 suggests that significant amount of NH3 observed in the blueshifted gas has 
been evaporated from dust grains. The evaporation of ammonia is supported by the fact 
that the ammonia abundance is enhanced in the blueshifted gas. If we assume that all of 
the NH3 molecules observed in the blueshifted gas were from dust grains, the observed ratio 



-10- 



of 1.3-1.7 suggest that the population of the rotational levels of NH3 at the time of ejection 
from the grain surface is represented by a Boltzmann distribution with the temperature of 
18-25 K. This formation temperature would be related to the excess energy distributed 
to the rotational degree of freedom in desorption processes, and are not necessarily equal 
to the dust temperature or the gas kinetic temperature (Takayanagi, Sakimoto, & Onda 
1987). It is most likely that the NH3 molecules retained on grains have been provided with 
sufficient energy to desorb from the grain surfaces by the passage of shocks (e.g. Nejad 
et al. 1990; Flower et al. 1995). The rotational temperature of 130-140 K suggests that 
the shock heating may be responsible for the desorption of the NH3 molecules. Sandford 
& AUamandola (1993) revealed that the sublimation of ammonia drastically increases as a 
function of temperature: the residence time of NH3 on a dust grain which is 10^^ yr at 40 
K becomes only 10~^ yr at 100 K. 

Recently, the infrared spectroscopic observations with Infrared Space Observatory 
(ISO) and ground-based telescopes reported that the abundance of the solid NH3 (relative 
to H2O) in icy grain mantle is no more than a few percent (e.g. Whittet et al. 1996a, b). 
These results imply that the NH3 molecule desorbed from grains is less important than 
previously expected (e.g. d'Hendecourt et al. 1985; Brown & Charnley 1990; Nejad et 
al. 1990). However, the ortho-to-para ratio measured in the L1157 outflow indicates that 
significant amount of NH3 arises from dust grains and that the gas-grain chemistry plays 
an important role in determining the NH3 abundance. 

We would like to thank the staff of NRO for the operation of our observations and their 
support in data reduction. We also thank Drs. S. Takano, K. Kawaguchi, Y. Hirahara, 
Y. Taniguchi, and J. Tkahashi for helpful comments. N.H. acknowledges support from a 
Grant-in-Aid from the Ministry of Education, Science, Sport, and Culture of Japan, No. 
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Fig. 1. — A map of CO {J— 1-0) outflow taken from Umemoto et al. (1992) with the observed 
positions marked. Three positions marked by flUed circles were observed with six transitions 
of NH3 from (1, 1) to (6, 6) and those marked by open circles were observed with the (1, 1) 
to (3, 3) transitions. Dashed contours indicate the distribution of the blueshifted emission 
(-5 to 1 km s~^) and solid contours show the redshifted emission (4.5 to 10 km s~^). 

Fig. 2.— NH3 spectra obtained at positions A (0", 0"), B (20", -60"), and C (45", -105"). The 
offsets are in arcseconds from the position of IRAS 20386+6751 [a(1950) = 20^ 38'"39."6, 
(^(1950) = 67°51'33"]. 

Fig. 3. — Rotation diagrams for the NH3 transitions measured toward positions B (20", 
-60"), and C (45", -105"). A^u; 5'u) ^-nd -Eu are the column density, the statistical weight, 
and the energy for the upper levels of the transitions. Filled circles denote the para species 
and flUed triangles are the ortho species. Vertical error bars give la errors. The solid lines 
between the (1, 1) and (2, 2) data points and the dashed lines between the (3, 3) and (6, 6) 
data points give the rotational temperatures derived from the combinations of the (1,1) and 
(2, 2) data and of the (3, 3) and (6, 6) data. If we assume an ortho-to-para ratio to be 1.7 
for position B and 1.3 for position C (the column densities of the ortho species are denoted 
by the open triangles), the data points of four transitions from (3, 3) to (6, 6) are fitted by 
the solid straight lines. 
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